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Abstract 

Two cases of a phenomenological model for ferromagnetism are considered, discrete and contin- 
uous. And the relationship, in general, between discrete and continuous models explored. In a 
similar way to the logistic map behavior, the continuous case is exactly solvable while the discrete 
one contains the bifurcation route to chaos. Through the ferromagnetic interpretation I comment 
on the relevance of this to understanding evolution of systems in time, the role of the configu- 
ration space in chaotic behavior, and how this understanding may lead to new exotic magnetic 
phenomena. 
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It has been known for some time in chaos theory that models behave in different ways 
when the dynamics is performed using different types of time evolution 0. A simple 
example of this is the logistic map, for which the continues time case is fully solvable for all 
finite parameter values while for discrete time the solution ranges from solvable to chaotic in a 
route known as the route to chaos. In this paper I analyze another simple example of this dual 
behavior, for a model based on two essential characteristics of ferromagnetism: symmetry 
breaking of the magnetization at the critical point and dissipating massless dynamics. The 
model considered is related to the Landau Ginzburg model of magnetism. By keeping the 
model simple the limit of the continuous case can be solved exactly and since the model 
contains the essential characteristics of ferromagnetism we have some tools that we can use 
to interpret the results. 

The outline of the article is as follows. The model is introduced and then analyzed with 
the assumption of continuous time evolution. For this case one finds an initial-condition 
dependent, predictable relaxation toward the fixed points of the potential energy, for all 
range of parameters. Next, the basis for discrete time evolution and the relationship between 
discrete and continuous time evolution are discussed. Finally the discrete case is explored and 
the same relaxation behavior is found for some parameter regime but new types of behavior 
(leading to chaos) for other parameter regimes. I finish by interpreting this difference, 
that results from a different treatment of time, and comment on the importance of this to 
our understanding of the way processes evolve with time, and indicate systems where such 
behavior may be found. 

The physical picture is as follows. Assume a bulk ferromagnet and let be its total 

free energy, with $ the total scalar magnetization in a given direction. If locally is a 

proper physical function, it must have an appropriate expansion as a power series. The sign 
of $ is an artifact of the choice of south magnetic pole, so = $) since reversing 

this choice should have no physical relevance. If we now take the fourth order approximation 
and neglect an additive term in the free energy we find, 

^=-V + ^$^ (1) 

with a, b the expansion coefficients. To ensure that infinite magnetization is excluded we 
must have 6 > 0. Since 6 is a second order parameter it is assumed for the rest that it is 
a constant. Equation (|l|) has two characteristic forms as we vary a, for a < there is a 
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single potential well around 0, we can identify this regime as paramagnetism. In the second 
regime, for a > 0, there is a double well potential. This broken symmetry regime is identified 
with ferromagnetism. Indeed for this model we can write a oc (Tc — T). Where T, Tc are 
the temperature and the Curie temperature respectively. 

Changing the magnetization is a strongly dissipative process that carries little to no 
inertia. To model this I use dissipative Newtonian dynamics for a massless particle and 
arrive at the following equation of motion for the magnetization 

6 = M(a$-6$^), (2) 

where M is the mobility and is related monotonically to a damping coefficient also 
known as the Khalatnikov damping constant. We can integrate this equation immediately 
to find the following solution for positive initial conditions, 
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For negative initial conditions we find the solution to be — <l>(t). For the trivial initial 
condition we find the trivial solution. The non-trivial long time behavior exhibits an asymp- 
tote located at the fixed point of the free-energy, limt_»oo ^(^) = R^{'^\J\)- This solution 
is completely deterministic for all finite values of parameters and initial conditions, by this 
it is meant that for any desired accuracy of prediction past some time to, one can give a 
requirement on the accuracy of the initial conditions at Iq which is good for all times t > tg. 
As we shall see shortly this situation is rather different in the discrete case. 

Before going into the discrete case we should explore the relationship between it and the 
continuous case. I claim there is a physical correspondence between the discrete and the 
continuous models. Usually the continuous solution would be called analytic and appear 
more pleasing to the trained physical eye. But there is no a-priori reason to assume that 
systems evolve continuously with time, passing through a continuous number of states. 
Indeed one can imagine how discrete processes are possible, where systems make jumps 
from state to state and between jumps remain in their state for some time interval At. 

To relate continues and discrete models one needs a mesh. For simplicity I assume a 
mesh with constant At, i.e. we have a discrete linear relationship for the magnetization. 
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= $(to + nAt). (4) 

In general an m order differential model will be given by the condition 
/($, <i>, . . . , <l>(™)(t)) = 0, with m initial conditions. Here, = A general rrfi^ 

order difference model is given by, gAti^n, ^n+i, ■ ■ ■ , ^n+m) = 0, with m initial conditions. 

Let gAt = be the difference model that one gets when discretizing a differential model 
/ = 0, by approximating the derivative ~ "^"t^^*^" similarly [^J for higher deriva- 
tives. And let C^t the discretizing transformation. Symbolically this can be writen as 
CAtif = 0) = {gAt = 0). From analysis we know that Cai — ^ J^d when At 0, where 
Id is the identity. From this general property one can show that for small enough At, the 
Cai mapping preserves the physical essence of the model. Stated more precisely, for a given 
mesh of time there is a unique correspondence between the differential and the difference 
models. Because we know that going from the difference model by taking the limit At — > 
one gets the corresponding differential model. If we show that this process is invertible, we 
can claim that the differential and difference equations of motion correspond to the same 
physical system with continuous or discrete evolution respectively. This invertiblity for 1** 
order systems is now shown. 

Assume for simplicity that a differential and a difference model are given by $ = 
/($), ^n+i = g{^n) respectively, where /, g are well behaved functions. For a given mesh 
size At > 0, let be given by equation (P and let CAt($ = /($)) = = fi'^n)) 

be the correspondence between differential and difference equations. Now if, 

CAt($ = /l($))=CAt($ = /2($)) 

it is easy to see that /i = /2, just by allowing the initial conditions to vary continuously, i.e. 
CAt is 1-1. To show onto, we can take /($) = -^[gi^) ~ ^] where g is the difference model. 
Then, 

CAt(^ = /($)) = ($n+l=^($n)). 

This shows that, for a given mesh of time, for every difference model there exists a unique 
differential model and vice versa. And thus, any changes that a given model will exhibit 
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when going from continuous to discrete time will result directly from this discretization and 
not from any other change to the underlying physical dynamics. 

I now show the discrete case for this model of magnetism. If we discretize equation (0) 
we find the following difference equation for the magnetization, 

= g{(^n) = $n + Ci[<l>„ - C2($n)^] (5) 

where Ci = Ma At and C2 = ^. At is a finite time element. 

The fixed points for this map, given by the condition $„+i = $„, , are $* = or 
$* = ±-FF- To understand the nature of these fixed points their stability is analyzed. 



^'($*) = l + ci[l-3c2(r)2]. 

Or, 

g'{0) = l + c,, 
g'{±J-) = l-2c,. 

V C2 

Thus for —2 < ci < the (j)* = fixed point is stable which corresponds to the para- 
magnetic state. For < Ci < 1, the $* = ^ stable fixed point for positive and 
negative initial conditions respectively. We recognize this as the ferromagnetic state. For 
larger values of Ci the fixes points become unstable and we find higher period fixed points. 
If we set Ci even larger we find a dense set of unstable periodic orbits with total measure 
zero This last scenario is known as chaos and the route we have taken is the period 
doubling route to chaos. 

Remembering that Ci = MaAt, M > and that for a forward evolving system At > 0, 
we can conclude again that a is the parameter that determines whether the system is a 
paramagnet or a ferro magnet. Note that for a backward evolving system (At < 0) the roles 
of paramagnetism and ferromagnetism change and we get the transition from ferromagnetism 
to paramagnetism as we lower the temperature. 

The standard conclusion derived from such a treatment is that as long as we keep At 
small enough such that ci is in the stable region (—2, 1) we can analyze equation with 
equation (Bl). But there is some thing more fundamental that we can conclude: the nature 



of time evolution - whether discrete or continuous, leads to different observation 
in physical systems. 

With this last conclusion in mind we are led quite directly to inferring from observations 
on the nature of ferromagnetism under discrete time evolution. For the system presented 
we consider two control parameters, the temperature and the mobility. Figure |I| is the 
bifurcation diagram as we vary the temperature, i.e. for a given temperature we iterate 
equation (H) many time steps (200) and we plot the final states it settles to (100 points). 
We start in the paramagnetic state with zero magnetization then at T^. (a = 0) we have 
a transition to the ferromagnetic state. As we keep increasing a we have more transitions 
until we reach the chaotic state. So, if time is discrete then we might expect that if we lower 
the temperature enough we will get a transition from the ferromagnetic state to a state 
for which the magnetization oscillates between two adjacent magnetization values. If we 
lower the temperature even further the magnetization difference between these two states 
will increase. However, it could turn out that time evolution is discrete and we would not 
observe such exotic magnetic states. The reason is that there is a physical bound on the 
value of a, set by the minimum possible temperature, T = Q K. It could turn out, especially 
if we expect a small At, that in most cases the corresponding value of Ci is well within the 
familiar ferromagnetic region. 

The second control parameter is the mobility. Figure ^ is the bifurcation diagram keeping 
a, h fixed, positive and varying the mobility. I remark that the part of the diagram for 
negative values of the mobility was left in mainly for esthetic reasons. As we increase the 
mobility we observe much the same bifurcation diagram as before. Indeed it is seen that 
there is a value of M for which the familiar ferromagnetic state ends and we start a new 
phase for which the magnetization oscillates between nearby states. Moreover there seems 
to be no natural cutoff for the mobility and detection of these exotic magnetic states seems 
more plausible Another system where detection of exotic states may be possible is 
discussed in 0. 

In the frame-work of magnetization these nonlinearities are intriguing and could also 
provide evidence about the nature of time evolution. Indeed as this and other work suggest, 
bifurcation routes and chaos rise naturally from discrete dynamics, and disappear once the 
continuous limit is taken. Interestingly enough, experimental and numerical work on a 
driven magnetic system do exhibit a period-doubling route to chaos [Q, ^ . Another possible 
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system where the exotic nature of the magnetization may be found is in a super-critical 
ferromagnet where the critical temperature (and thus a) is very large compared to everyday 
ferromagnets. In such systems the distance between nearest neighbor magnetic moments 
is small such that the overlap between adjacent wave functions and hence the exchange 
integral are greatly increased. An example of such a system is a neutron star, where the 
distances between adjacent magnetic moments are of nuclear order. Such stars exhibit a 
large-scale ferromagnetic ordering at very high temperatures which gives a very high critical 
temperature. When such a star cools down a becomes large and the possibility of exotic 
states arises in the realm of discrete time. Such exotic states will oscillate between different 
magnetization states with period of the order At. If our instruments make observations over 
a time much larger than At we would observe only the average magnetization. 

In conclusion, a simple model for the bulk magnetic state has been studied. It is found 
that under continuous time evolution assumption this model admits only a steady-state an- 
alytical solution that is fully deterministic for all model-parameters values. Under discrete 
time evolution assumption the model exhibits the same analytical solution for some parame- 
ter values but as part of the period doubling route to chaos. This different behavior resulted 
from the different assumption with regard to time evolution and thus such observation can 
aid in understanding the nature of time and of time evolution. 

The author would like to acknowledge the support generously provided by Natali J. 
Teszler. 
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FIG. 1: Temperature bifurcation diagram with constant mobility for the total magnetization in the 
discretized Ginzburg-Landau model, MAt = 1 and 6=1. Top: Positive initial condition Bottom: 
Negative initial condition. 



9 




FIG. 2: Mobility bifurcation diagram with constant temperature for the total magnetization in the 
discretized Ginzburg-Landau model, a At = 1 and C2 = 1. Top: Positive initial condition. Bottom: 
Negative initial condition. 



10 



